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Transfer-matrix methods are used for a tight-binding description of electron transport in graphene¬ 
like geometries, in the presence of spin-orbit couplings. Application of finite-size scaling and phe¬ 
nomenological renormalization techniques shows that, for strong enough spin-orbit interactions and 
increasing on-site disorder, this system undergoes a metal-insulator transition characterized by the 
exponents v = 2.71(8), r] — 0.174(2). We show how one can extract information regarding spin 
polarization decay with distance from an injection edge, from the evolution of wave-function am¬ 
plitudes in the transfer-matrix approach. For (relatively weak) spin-orbit coupling intensity fi, we 
obtain that the characteristic length As for spin-polarization decay behaves as As oc 
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I. INTRODUCTION 


results; finally, concluding remarks are made. 


In this paper we consider electronic transport in two- 
dimensional (2D) honeycombe-lattice structures with 
spin-orbit (SO) interactions. While single-parameter 
scaling predicts that noninteracting electron states in 
2D are generally localized by disorder—, it is known 
that in this marginal dimension the enhancement of for¬ 
ward scattering provided by SO effects can partially off¬ 
set disorder-induced quantum interference, thus giving 
rise to a conducting phase^Ti^. Motivation for consid¬ 
ering a honeycomb geometry is provided mostly by re¬ 
cent progress in experimental synthesis, and theoretical 
understanding, of graphene^^— . It should be noted that 
some of carbon’s cousin elements in group IVB have been 
predicted and, to varying degrees of success, shown to 
crystallize in a similar honeycomb arrangement, giving 
rise respectively to silicene^ and germanene^. Stable hon¬ 
eycomb structures have also been found for boron nitride 
(h-BN)i£ as well as Si02 (silicatene)^^ and other elements 
or compounds. The ideas and methods exhibited in what 
follows are, in principle, applicable to any of the above. 
However, where comparison of our results to experimen¬ 
tal data is pertinent we generally restrict ourselves to 
graphene, as this is so far the best-understood member 
of the ensemble. 

Our purpose here is twofold: (i) to investigate location 
and universality properties of the second-order metal- 
insulator found in these systems, for strong SO coupling, 
with increasing on-site disorder; and (ii) to examine the 
spatial evolution of the polarization of an initially fully 
spin-polarized current, injected into one such hexagonal- 
lattice system with SO couplings. 

Sectionllllbelow recalls selected existing results, as well 
as some technical aspects of the transfer-matrix (TM) 
method used in our calculations. In Sec. m we discuss 
the metal-insulator transition. In Sec. IIVI we show how 
one can extract information on the evolving state of po¬ 
larization of an injected electron beam, from the anal¬ 
ysis of wavefunction amplitudes generated via the TM 
method. In Sec. El we provide a global analysis of our 


II. THEORY 


The model one-electron Hamiltonian for this problem 
can be written as 
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where Cia are creation and annihilation operators 
for a particle with spin eigenvalues tr = ±1 at site i, 
and the self-energies Si are, in general, independently- 
distributed random variables; Vijaa' denotes the 2x2 
spin-dependent hopping matrix between pairs of nearest- 
neighbor sites whose elements must be consistent 

with the symplectic symmetry of SO interactions. 

Although modeling SO coupling in real systems usu¬ 
ally relies on including e.g. p orbitals, with their corre¬ 
sponding degeneracies, in Eq. (HD we resort to the cus¬ 
tomary description via an effective Hamiltonian with a 
single (s—like) orbital per site; for a pedagogically clear 
explanation of how this approach works from first princi¬ 
ples, see e.g. the Appendix of Ref. 0. Since we shall not 
attempt detailed numerical comparisons to experimental 
data, such formulation seems adequate for our purposes. 

In the effective-Hamiltonian context, there is some lee¬ 
way as the specific form of the hopping term is con¬ 
cerned, depending on whether one is specifically con¬ 
sidering Rashba- or Dresselhaus- like coupling s ^ i - — , or 
the focus is simply on the basic properties of systems 
in the symplectic universality classi^^— . One constraint 
is that it must incorporate the basic symmetries of SO 
coupling, which have close connections with quaternion 
algebrapiSiii. 

Here we use the implementation of Refs, [l^ and 
namely: 
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where / is the 2x2 identity matrix, cr* are the Pauli 
matrices, and fj, gives the intensity of the SO coupling; 
below we consider the (real) {P^} as either uniform 
(V^ = yy = on all bonds) or randomly (uniformly) 
distributed in [—1/2,1/2]. Thus all energies will be writ¬ 
ten in units of the /i = 0 nearest-neighbor hopping. 

Similarly the £i will either be constant, site- 
independent, or taken from a random uniform distribu¬ 
tion in [—W/2,W/2\. 

The form Eq. ([5]) for the hopping term does not exhibit 
the explicit multiplicative coupling between momentum 
and spin degrees of freedom, characteristic of Rashba-like 
Hamiltonians^ii^. A similarly-decoupled effective Hamil¬ 
tonian can be derived from first principles for carbon 
nanotubes, see Eqs. (3.15), (3.16) of Ref. [l^ In two 
dimensions one should not expect significant discrepan¬ 
cies between results from either type of approach, as long 
as one is treating systems without lateral confinement. 

We apply the TM approach specific to tight-binding 
Hamiltonians like Eq. to finite-width strips of the 

honeycomb lattice, with N sites across. Adaptation from 
the more usual square-lattice geometry is straightfor¬ 
ward, closely following the lines used for the TM descrip¬ 
tion of localized (e.g., Ising and Potts) spin systems22ri24. 
Two distinct orientations are possibl o^^'^^ , with the TM 
proceeding either (a) perpendicularly^^ or (b) parallel to 
one lattice direction^^. Case (a) corresponds to a "brick" 
lattice, i.e., a square lattice with vertical bonds alter¬ 
nately missing. 

In the terminology of quasi one-dimensional carbon 
nanotubes (CNT) and nanoribbons (CNR)^, a narrow 
strip with periodic boundary conditions across in geome¬ 
try (a) would be topologically equivalent to an armchair 
CNT, while one in geometry (b) would correspond to a 
zigzag CNT. Conversely, free boundary conditions paral¬ 
lel to the TM’s direction of advance give: zigzag CNR in 
(a), armchair CNR in (b). 

Detailed consideration shows that implementation of 
the TM scheme of Ref. [l^ in geometry (b) involves a 
number of cumbersome intermediate operations [ mostly 
matrix inversions, see Eqs. (8)-(16) of Ref. [2^]. In what 
follows we always make use of geometry (a), for simplic¬ 
ity. 

We briefly recall selected aspects of the TM formula¬ 
tion introduced in Ref. [11 and of its adaptation for a 
honeycomb geometry with SO couplings. 

Consider a strip of the square lattice, cut along one 
of the coordinate directions. For the orthogonal univer¬ 
sality class with site disorder, denoting by A: = 1,..., M 
the successive columns, and i = 1,...7V the respective 
positions of sites within each column of a strip, the re¬ 
cursion relation for an electronic wave function at energy 
E is given in terms of its local amplitudes (which can all 
be assumed real), {aife(E)}, and tight-binding orbitals 


\ik), as: 
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where / is the N x N identity matrix, the energy depen¬ 
dence has been omitted for clarity, and (invoking periodic 
boundary conditions across the strip), 

/E — Elk —1 0 • • • —1 

— “I E — £2k — 1 ••• 0 

\ —1 0 ... —1 E — £]sik 



For a honeycomb lattice in the "brick" geometry, the 
changes to Pk (i) the off-diagonal elements are 

of the form — (l -|- (—1)*"'’^)) /2, reflecting the alternately 
missing vertical bonds; and (ii) an elementary step con¬ 
sists of applying the TM twice, in order to restore peri¬ 
odicity. 

The introduction of SO couplings along the bonds, see 
Eqs. ([T|), (U), means that the Qik are now spinors, written 
on the basis of the eigenvectors of cr^ as: 

(5) 

where the a/), are complex. The matrices Pk and 
the subdiagonal I of Eq. ([3]) are now 2N x 2N, while the 
diagonal terms of Pk are doubly degenerate, and the non¬ 
zero off-diagonal ones are replaced either by the (bond- 
dependent) negative of matrix Vij of Eq. ([3]) or that of 
its hermitean adjoint , depending on whether they are 
supra- or sub-diagonal21. The negative unitary elements 
of the (now 2N x 2N) supradiagonal identity matrix of 
Eq. ([3]) are replaced by the bond-dependent negative of 



III. METAL-INSULATOR TRANSITION 
A. Introduction 

To make contact with previous work on the square 
latticei^^— , we initially considered systems without site 
randomness, i.e., all £i = 0, and two versions of SO cou¬ 
pling: (i) uniform with /r = 0, = 1 

in Eq. o, ( ii) random, with = 2 (as in Ref. El) and 
R^, R*^, R^ uniformly distributed in [—1/2,1/2]. 

In case (i), analysis of the resulting dispersion re¬ 
lation shows that the allowed states occupy a band 
with the same structure as that for a system with 
/r = 0 but, analogously to the square-lattice sys¬ 

tems with uniform SO term of Ref. Il4l . in the range 
±3^/1 + [(R-)2 + (Ry)2 + 

For case (ii) the corresponding density of electronic 
states (DOS) p(e) can be evaluated by making use of 
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eigenvalue-counting theorems^^— . Our implementation 
takes advantage of the sparse nature of the Hamiltonian 
matrix written on the site basis, and resorts to Gaussian 
elimination algorithms on strip-like geometries (with pe¬ 
riodic boundary conditions across), closely following the 
steps described in Refs. [3ll - l^ . 

Since p(e) is calculated from the finite difference be¬ 
tween successive values of integrated DOS up to adja¬ 
cent energies e and e + A2i“— , the bin size A has to be 
optimized in order to reduce oscillations in the numer¬ 
ical results while still capturing relevant structural de¬ 
tails of the DOS. We generally took strip-like systems 
with N > 80 sites across, and length M > 200, for 
which A = 0.06 proved to be a reasonable choice. For 
cases such as (ii) where quenched randomness can play 
a role in inducing further fluctuations, we saw that for 
the system sizes used, the self-averaging provided by hav¬ 
ing a large(ish) number of local disorder realizations was 
enough to render such effects relatively unimportant. 

We numerically evaluated p{e), for both (i) [as an in¬ 
dependent check of the soundness of our algorithms ] and 
(ii). Results are shown in Fig. [T] For (i) we used po = 
l/-\/3, which ensures that + (F^)^ + (F^)^] 

equals p? [((F"^)^ -I- (F^)^ -I- (F^)^)] of case (ii), where 
angular brackets stand for ensemble average. 

One sees that the resulting bands indeed have the same 
width, although the shape of tails at the edges differs. As 
expected, these are rather abrupt in (i) and smoother in 
(ii), in line with the probabilistic character of the latter’s 
coupling distribution. Furthermore, while the band in 
case (i) keeps all qualitative features of the p = 0 system, 
including the zero at e = 0 , these are lost in case (ii); for 
example, the shallow minimum at the origin corresponds 
to p(0) = 0.077(2). 

Note that the effective strength of intrinsic SO cou¬ 
pling in graphene is estimated to be 25 — 50 peV—, of 
order 10“^ of the nearest-neighbor hopping 70 = 2.8 eV— . 
The values used in this Section are for illustration only, 
and not intended to reflect the actual properties of pure 
graphene. We return to this point in Sec. IVl below. 

B. Phenomenological renormalization 

Here the SO coupling is represented as in model (ii) 
of Sec. nil A1 with p = 2 (hxed). We now introduce site 
randomness, i.e. P(ei) = constant, |ei| < W/2 and zero 
otherwise, and allow the respective distribution width W 
to vary. 

Following standard proceduresi^ we considered E = 0, 
at the Dirac point of the unperturbed Hamiltonian, and 
iterated the TM on strips of width N sites and length 
M ^ N with periodic boundary conditions across. The 
characteristic Lyapunov exponents were extracted22i2i, 
with the longest localization length Xn being given by 
the inverse of the smallest of those. According to finite- 
size scaling and the phenomenological renormalization 
ansatz, we plotted the scaled localization lengths A^r = 
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Figure 1. (Color online) Density of states (DOS) for the 
Hamiltonian Eq. Q. Dashed line, red: p = 0, exact re¬ 
sult. Curves (i) and (ii) evaluated by eigenvalue-counting on 
300 X 300 systems. Curve (i), magenta: uniform SO coupling, 
po = l/%/3. Curve (ii), blue: random SO coupling, p — 2 (see 
text). 

Aat/A against varying W, looking for the mutual in¬ 
tersection of the curves corresponding to different strip 
widths which gives the location of the critical point of 
the metal-insulator transition. 

We took N = 10, 12, 14, 16, 24, 32, and 48, with 
M = 10®. The results are shown in Fig. 

It is known that nonlinearity of scaling fields and/or 
irrelevant variables^ can cause sizable distortions in the 
estimates of critical quantities at Anderson localization 
transitions. Efficient procedures have been devised to 
correct for effects of this sort 15,36,37 ^ noted that if 
only strip widths A > 10 were used, corrections to scaling 
due to irrelevant fields (usually dealt with by methods 
explained in Refs. and [ 3 ^ were of little import. Thus 
we concentrated on accounting for nonhnearitiesi^. 

Following the lines of Ref. [H [see especially their 
Eqs. (6)-(ll)[, we dehne Wc as the disorder distribu¬ 
tion width at criticality, and with w = (IF — IFc)/IFc, we 
assume that: 

InAjv = , (6) 

where ip is a smooth function of E and IF which goes 
to zero at criticality. At fixed E, one then considers a 
truncated Taylor expansion of ip in terms of w: 

Tlfp 

ip = '^ipkw^ . (7) 
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Figure 2. (Color online) Systems with random SO coupling, 
p = 2: raw data for scaled localization lengths against width 
of site-disorder distribution. Uncertainties are of order of 
symbol sizes or smaller. 


Figure 3. (Color online) Systems with random SO coupling, 
p = 2: scaling plot of (logarithm of) scaled localization 
lengths against |^/)| [see Eqs. 1- Key to symbols is the 

same as in Fig. [2] 


Table I. Adjusted values of parameters for the scaling plot of 
Fig. (3] with 'i/ji = 1 (fixed) in Eq. Q. 


p 

2.71(8) 

In Ac 

0.459(9) 

Wc 

6.21(2) 

1p2 

-0.06(1) 

ai 

-0.158(5) 

a2 

-0.0076(12) 

as 

0.00028(8) 


Nonlinearities in the argument of F are thus taken into 
account with > 1. Plugging this back into Eq. ([5]) 
one gets another truncated Taylor series near criticality: 

mo 

F{x) = In Ac + ^ amx"^ , ( 8 ) 

m—1 


accurate estimates, mostly in the range 2.6 — 2.9 [see 
Table 1 of Ref. 1^ ]. 

For comparison of In Ac = 0.459(9) with results per¬ 
taining to the square lattice, one must recall the geo¬ 
metric correction factors of the honeycomb lattice 22 r^. 
This means that, in order to produce an estimate of 
the decay-of-correlations exponent 77 = (ttAc)”^ given 
by conformal invariance^^, the raw TM result for Ac 
must be multiplied by a factor of 2fy/3. We thus ob¬ 
tain = 1.826(15) [so ?7 = 0.174(2) [which com¬ 

pares rather well with Ac = 1.844(4) of Ref. We note 
that the critical-amplitude results of Ref. do not com¬ 
pare directly with ours because those authors study the 
scaling behavior of the typical localization length. As is 
known, this quantity is given by the zeroth moment of the 
correlation-function probability distribution^, whereas 
here we deal with average quantities, i.e. ones related to 
the corresponding first-order moment. 


where Ac = Ajv-).oo(7?, w = 0). One can either set ipi = ^ 
in Eq. ©, or oi = 1 in Eq. ([5]) without loss of generality. 
It is expected that by using the logarithm in Eq. ®, the 
smoothness assumption underlying the Taylor expansions 
will be fulfilled with fewer terms than if the Ajv them¬ 
selves were considered. 

Fig. [3] shows the best-fitting scaling plot of the data of 
Fig. for which we took njjj = 2, uiq = 3. Numerical 
values are given in Table HI 

The adjusted v — 2.71(8) is in very good agreement 
with v = 2.75(4) of Ref. which is, to our knowledge, 
the most accurate result to date . Earlier work gave less 


IV. SPIN RELAXATION 

In studies of spintronics in semiconductor a^^d^ the 
spin coherence length is one of the quantities of inter¬ 
est. Here we consider the decay of spin polarization in 
electronic transport in quasi one-dimensional geometries, 
in the presence of SO couplings. The subject is usu¬ 
ally approached via Green’s function techniquesi^. We 
show that this problem can be investigated by consider¬ 
ing the evolution of wavefunction amplitudes in the TM 

contex t 
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Usual TM treatments focus on extracting the spectrum 
of characteristic Lyapunov exponents, which demands re¬ 
peated iteration along M ^ 1 columns, with frequent 
orthogonalization to avoid contaminatiooi^Ti^i. However, 
one sees that the recursion relation synthesized in Eq. m 
contains information on how the electron wavefunction 
evolves, starting from specified initial conditions. In pres¬ 
ence of SO couplings the off-diagonal hopping matrix el¬ 
ements induce spin flips, thus affecting spin polarization 
along the system’s length. 

Assume that a fully spin-polarized electron beam is 
injected into the system, i.e., uJq, ^ 0, a|g = = 0, 

i = One can extract information about the 

beam’s polarization state M lattice spacings down the 
strip by iterating Eq. m M times and examining the 
resulting coefficients {a^^}. In this case the beam 

polarization V at column M is given by: 


V{M) 




(9) 


The initial conditions just mentioned can be viewed 
as representing an ideal lead, i.e., one without SO inter¬ 
action, from which the beam is injected into the strip. 
Although, for a complete description one would need to 
take into account reflection at the injecting boundary (as 
well as at the strip’s end, presumably linked to a second 
ideal lead), these features do not influence the calculated 
polarization decay length, as this quantity is computed 
from ratios of (sums of squared) amplitudes each taken 
at a fixed positioi>i^. 

We evaluate spin polarization profiles on systems with 
periodic boundary conditions across, i.e., topologically 
equivalent to CNTs; we keep E = 0, fixed, as this corre¬ 
sponds to the Eermi energy which is the relevant level for 
transport phenomena. The SO couplings are again rep¬ 
resented by model (ii) of Sec. IIIIA1 although the overall 
amplitude ^ will be allowed to vary. Our results are aver¬ 
ages over typically 10^ independent samples. For each of 
these we generate random sets of the {aJoi^ 
as the {Vij}. Thus we are sampling over the ensemble of 
steady-flow configurations. 

We consider the simplest case with no site random¬ 
ness, i.e., all ETj = 0 in Eq. ©■ This removes one (finite) 
length scale from the problem, as one would have the 
[site] disorder-associated mean free path Ay —>■ oo. 

Fig.[4]shows that spin polarization P(x) settles into ex¬ 
ponential decay against position x along the nanotube’s 
axis, after a short transient region of steeper variation. 
The characteristic decay lengths A^ are found by adjust¬ 
ing the appropriate sections of numerical data to a pure 
exponential-decay form. 

For fixed /i, we found a slight dependence of As on N, 
which is strongest for small N < 20. As N increases, 
saturation becomes evident and it is possible to estimate 
As (/r) = limjv-j.oo As (p., N) with good accuracy by using 
data for N up to 100 (see the inset to Fig. B- 

The main diagram in Fig. [5] shows that, to a very good 


X 





Figure 4. (Color online) Log-linear plot of spin polarization 
Vix) versus position x along axis of nanotube with A = 30 
sites across. A fully-polarized beam is injected at a: = 0 with 
E = 0. Systems with random SO coupling strength p as 
shown, site-disorder width W — 0. The characteristic decay 
lengths As are adjusted values from fits of data for a: > 30 
(fitted lines are shown in color). 


extent, the relationship As(/r) oc holds for the range 
0.05 < fi < 0.25. This is in line with the findings of 
Ref. [13 for quantum wires with nonzero on-site disorder 
and Rashba-like SO coupling. Such inverse-square depen¬ 
dence of As on /i is thus likely to be a universal property 
for electronic transport in two-dimensional systems with 
SO effects. 

Next, we introduce quenched impurities with SO inter¬ 
action onto an otherwise pure system, i.e. one in which 
such interaction is generally absent. This amounts to a 
simplified representation of pure graphene (where, as re¬ 
called in Sec. IIII Al SO coupling is very weak) doped with 
suitable impurities. 

Many experimental and theoretical studies deal with 
impurities on CNRs, where edge effects play an impor¬ 
tant role in the energetics of favored defect locations . 
By considering only nanotube geometries, here we need 
not account for this sort of positional preference inhomo¬ 
geneity. Furthermore, we restrict ourselves to so-called 
hole defects, i.e. adatoms which sit over the center of a 
hexagon^. 

Denoting by p the fraction of randomly chosen 
hexagons which have an impurity atop their center, we 
assume that hopping along all six bonds which make such 
a hexagon is characterized by SO couplings Vij as given 
in Eq. (21), with the V^, V^, randomly distributed as 
in case (ii) of Sec. IIII Al We take /x = 2 as in Sec. IIIII 
which makes the adatom’s average SO interaction about 
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Figure 5. (Color online) Main diagram: double-logarithmic 
plot of As(/r) = limjv-^oo As(/r, against random SO cou¬ 
pling strength fi, for site-disorder width W = 0 and energy 
E = 0. Full line is As(/i) = 3.0 (best fit to data). Inset: 

As{^J.,N) against 1/A^ for ^ = 0.11. 

as strong as the hopping term of pristine graphene. 

We consider only p < 0.01 and neglect effects due to 
adjacent impurities, which should be a reasonable ap¬ 
proximation in this concentration range. 

In Fig.iniit can be seen that for p = 0.005 the relaxation 
length As is very close to that for a system with SO 
coupling on all bonds, and p, = 0.167. For relatively 
short distances x <20 from the injection edge the initial 
decay rate is found to be distinctly higher for the latter 
case than for the former, as highlighted in the inset of 
Fig. We refrain from ascribing much significance to 
this difference since we do not expect our approach to 
give an accurate description of such short-range effects. 
So, concerning the region farther than some 20 lattice 
spacings from the origin, the effect on (asymptotic) spin 
polarization decay of 3% of bonds with SO coupling p = 2 
is similar to that of all bonds having p — 0.167. 

We checked the dependence of A^ on p and p. For fixed 
p = 2 we took p = 0.003, 0.005, 0.0075, and 0.01 which 
gave a rather good fit to a p~^ dependence with x = 
1.1(1); then for fixed p = 0.0075 we additionally made 
p = 0.75, 1.0, and 1.5. In this case one gets As oc 
with y Ri 1.7. If one instead assumes that y = 2, found 
for smaller p < 0.25 and systems with SO couplings on 
all bonds (see Fig. 0, still applies here, one gets: 

Asip,p) = ^, A = 2.4(5), (10) 

PfJ- 

where the large scatter in the estimate of A reflects 
the aforementioned poor quality of the fits of behavior 


Figure 6. (Color online) Main diagram: log-linear plot of spin 
polarization Vix) versus position x along axis of nanotube 
with = 30 sites across. Hexagons (red): SO coupling p = 2 
on a fraction p = 0.005 of hexagons, zero elsewhere (see text); 
squares (blue): SO couplng p = 0.167 on all bonds. Values of 
As from fits of a: > 40 data to exponential form. Inset: same 
data, linear scale on vertical axis, close-up of a; < 20 region. 

against varying p. Nevertheless, Eq. m can be a rough- 
and-ready guide to estimate the relaxation length for val¬ 
ues of p and p closer to physically realizable experiments 
than those used here. 


V. DISCUSSION AND CONCLUSIONS 

We have studied model tight-binding Hamiltonians for 
the description of electron transport in graphene-like ge¬ 
ometries, in the presence of SO couplings. Our main 
interest is in the behavior of very large sheet-like sam¬ 
ples, although in the quasi-one dimensional context of the 
TM methods applied here, the use of periodic boundary 
conditions across makes our systems topologically akin 
to CNTs. For the latter type of system, especially in 
the limit of very narrow nanotubes, a more realistic de¬ 
scription should probably incorporate curvature effects, 
as done in Ref. HH 

The energy unit used here equals the nearest-neighbor 
hopping parameter in pristine graphene, thus translating 
into 2.8 eV— . The intrinsic SO coupling in graphene is 
estimated as 25—50 peV— ; weakly-hydrogenated samples 
have been reported as giving a colossal enhancement on 
this by three orders of magnitude^. The latter value 
would then correspond to p r; 10 “^, similar to the low 
end of the range investigated in Sec. lIVl for systems with 
SO couplings on all bonds. On the other hand, we used 
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/j. = 2, both in Sec. IIIII and for impurity adatoms in 
Sec. lIVl partly for comparison with extant worbi^^— , and 
also in order to produce well-defined numerical results 
amenable to unequivocal interpretation. 

In Sec. ini we showed that for strong SO interactions 
a second-order metal-insulator transition takes place, 
which is in the same universality class as that found 
for square-lattice systems^d^^— . Note that the energy 
E = 0 used in our calculations corresponds to the Dirac 
point of pure graphene, at which (contrary to the square- 
lattice case) the tight-binding DOS vanishes identically. 
Although in real graphene a (quantized) non-zero mini¬ 
mum conductivity has been found at the Dirac point^, 
the effect seen in our results is unrelated to this, being 
of a larger order of magnitude (the metallic phase ex¬ 
tends up to site disorder of strength Wc ~ 6.2). In fact, 
as can be seen in Figure [U the random SO couplings 
used in the model investigated in Sec. IIII Bl account for 
the significant departure from zero of the DOS close to 
E = 0. So, it is the latter feature which sets the stage 
for the relative robustness of the conducting phase in our 
model. The above remark should apply also to the case 
of structures mentioned in Sec. HI such as boron nitride^ 
or silicateneii, which are wide gap insulators. Indeed, 
it is expected^ that in general the introduction of ran¬ 
domness will be accompanied by the appearance of states 
outside the pure-system bands, in particular within the 
gap. Whether that will be enough to give rise to a con¬ 
ducting phase should depend on quantitative details of 
the disordered potential. 


In Sec. IIVI we showed how one can extract information 
regarding spin polarization decay with distance from an 
injection edge, from the evolution of wave-function am¬ 
plitudes in a TM context. We illustrated the pertinent 
ideas in the simple context of a nanotube-like geometry, 
with no on-site disorder, and investigated the dependence 
of the spin relaxation length Ag on SO coupling strength 
/i. For small p, closer to physically realizable values, 
we found the dependence Ag oc which seems to be a 
universal relationship for two- (or quasi-one) dimensional 
systems with SO interactionsi^. For SO couplings acting 
only on impurity sites, randomly distributed with low 
concentration p ^ 1, we found numerically Ag cx p~^, in 
line with elementary probabilistic considerations. 

It must be noted that modeling SO coupling via 
Eq. Q with the V^, k = x,y,z randomly distributed in 
[—1/2,1/2] gives an effect which is, on average, isotropic 
in spin space. So this formulation is not suitable, e.g., 
for a realistic discussion of precession effects. 

Prospects for further application of the ideas presented 
in Sec. lIVI would include taking site disorder into account, 
as well as studying CNR geometries. For the latter, in¬ 
homogeneities in local current density in the transverse 
direction to average flow would be directly accessible. 
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